Direct photons from relativistic heavy ion collisions at CERN SPS and at RHIC 



(N 
O 

o 

(N 



> 

o 
m 
o 

O 



o 

=5 



X 



A. K. Chaudhurifcl 

Department of Physics, Technion-Israel Institute of Technology, 
Haifa 32000, Israel 
and 

Variable Energy Cyclotron Centre 
l/AF,Bidhan Nagar, Calcutta - 700 064 

T. Kodama 

Instituto de Fisica, Universidade Federal do Rio de Janeiro 
Cxp 68528, Rto de Janeiro, 21945-970, RJ, Brazil 



Assuming QGP as the initial state, we have analyzed the 
direct photon data, obtained by the WA98 collaboration, in 
158 A GeV Pb+Pb collisions at CERN SPS. It was shown, 
that for small thermalisation time, two loop rate contribute 
substantially to high pr photons. We argue that for extremely 
short thermalisation time scale, the higher loop contribution 
should not be neglected. For thermalisation time 0.4 fm or 
greater, when higher loop contribution are not substantial, 
the initial temperature of the QGP is not large and the system 
does not produce enough hard pr photons to fit the WA98 
experiment. For initial time in the ranges of 0.4-1.0 fm, WA98 
data could be fitted only if the fluid has initial radial velocity 
in the range of 0.3-0. 5c. The model was applied to predict 
photon spectrum at RHIC energy. 

PACS numbers(s):12.38.Mh,13.85.Qk,24.85.+p,25.75.-q 



Ever since possible existence of deconfining phase tran- 
sition from hadrons to quark gluon plasma (QGP) was 
predicted by the Quantum Chromodynamics (QCD), nu- 
clear physicists are eagerly looking for the experimen- 
tal evidences of QGP, the new state of matter |EJ. It 
has even been claimed that, evidence of the new state 
of matter is already seen at CERN SPS in Pb+Pb col- 
lision The claim was based mainly on the experi- 
mental J /if) suppression and strangeness enhancement, 
in Pb+Pb collisions. Due to intrinsic ambiguities in dis- 
entangling the similar effects arising from hadronic pro- 
cesses, their claim is rather considered controversial ||. 
On the other hand, it may be noticed that electromag- 
netic signals (dileptons and photons), were not consid- 
ered in establishing the claim, even though, photons and 
dileptons are well recognized probes of QGP. They suffer 
minimal final state interactions and can give information 
about the early stage of the collisions. Experimentally 
however, it is a challenge to obtain precision photon data 
from huge background. Recently WA98 collaboration has 
published their single photon emission data for 158 A 
GeV Pb+Pb collisions at CERN SPS {§. Several au- 
thors have analyzed the WA98 single photon data , 
but their conclusions are not convergent. It seems that 
the data could be explained in a QGP or in a hadronic 
gas scenario, with or without initial radial velocity. 



In the present paper, we further analyze the WA98 
data in the initial QGP scenario. We first assume that 
the thermalisation time Tj, beyond which hydrodynamics 
is applicable, is larger than 0.4 fm. The reason for this 
is the uncertainty in the elementary photon production 
rate. For smaller tj, higher loops contributions can be 
substantial, ruining any agreement obtained with data. 
If Tj > 0.4 fm is the case, we then confirm the analysis 
of Peressounko et al || that the data require an initial 
radial velocity. It was also shown that the high pt pho- 
ton spectra depend quite strongly on the initial radial 
velocity distribution. For realistic (surface peaked) ve- 
locity distribution, initial radial velocity in the range of 
0.3-0. 4c is required. 

Procedure for obtaining photon spectra in a hydro- 
dynamic evolution is well known |lfj . We have solved 
the hydrodynamic equations <9 A1 T A " / = for a baryon free 
gas assuming cylindrical symmetry and boost-invariance. 
The equation of state for QGP was assumed to be p q = 
a q T 4 — B with a q — 42.257r 2 /90. The hadronic equation 
of state was generalized to include all the mesonic reso- 
nances with mass < 2 GeV. The cut off 2 GeV is rather 
arbitrary and we verify that the results do not depend on 
the value of cut-off significantly. The bag constant B was 
obtained from the Gibbs condition pqgp(T c ) = Phad(T c ). 
The critical temperature (T c ) and the freeze-out temper- 
ature (T F ) were assumed to be 180 MeV and 100 MeV 
respectively. 

Hydrodynamic models require, as inputs, the initial 
time (Tj) and initial energy density or equivalently the 
temperature (Tj) of the fluid. Tj is the thermalisation 
time, beyond which quasi-equilibrium is established and 
hydrodynamics is applicable. Tj and Tj are parameters 
of the model, unless one obtains them from some micro- 
scopic transport models. In order to reduce parameters, 
one further assume that the fluid flow is isentropic. Then 
for a given Tj the initial temperature Tj of the fluid (QGP) 
can be obtained by relating the entropy density with the 
observed pion multiplicity (assuming pion decoupling to 
be adiabatic) JTl| l, 

where c = 27r 4 /45£(3) and Ra is the transverse radius 



of the system (assumed to be 6.4 fm for Pb+Pb colli- 
sions). 6 = corresponds to central collisions. In table 
1, initial temperature of the QGP, for certain values of 
the initial time Tj are shown. Rapidity density was as- 
sumed to be dN/dY— 750 ||. It can be seen that choice 
of initial time is of crucial importance. QGP is formed at 
higher temperature for low values of the thermalisation 
time. High px photon production rate depends strongly 
on the initial temperature and it is possible to increase 
its yield by arbitrarily reducing the thermalisation time 
and correspondingly increasing the initial temperature of 
the QGP fluid. 

For the single photons from hadronic gas we include 
the following processes, 

(a) mr — > p7, (b) irp — > irj, (c) u> — > irj, (d) p — > ttttj 
(e) irp — ► A\ — > 7T7 

rates for which are well known [T^jr^ |. 

Rate of production of hard photons from QGP were 
evaluated by Kapusta et al %M- To one loop order, 
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where £ ~ 0.23 is a constant. The summation runs over 
the flavors of the quarks and e/ is the electric charge of 
the quarks in units of charge of the electron. 

Recently Aurenche et al evaluated the production 
of photons in a QGP. At two loops level Bremsstrahlung 
photons {qq{g) — ► qq{g)~1 found to be dominating the 
compton and annihilation photons. 

The rate 
of production of photons due to Brcmsstrahlung (cor- 
rected for the factor of 4 |l^,|l6|) was evaluated by them 



as, 
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E^- = -^aa s > ' e 2 f —e~ E/T {J T - J£)I{E,T) (3) 



where Jt ~ 4.45 and Jl ~ —4.26 for two flavors and 
3 colors of quarks. For 3 flavor quarks, Jt ~ 4.8 and 
J L 4.52. I{E,T) stands for, 

I(E, T) = [3C(3) + ^| + (|) 2 In 2 + AU^-e^) 

+2Li 2 (-e-W T ) - {E/Tf ln(l + e -^ T )} (4) 
and the poly-logarith functions Li are given by, 
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Aurenche et al jl5| also calculated the contribution of 
the qq with scattering, which should also be corrected for 
the same factor of 4. The corrected rate is, 

E ^ p = H aa ^ ET ^ E/T ( J T-J L ) (6) 



We would like to call attention to the fact that two- 
loop photon rate from QGP is not complete. Higher loops 
contribute to the same order pl|. Also Landau-Migdal- 
Pomeranchuk effect has been neglected Questions 
naturally arise how can then one obtain a reliable es- 
timate of photon yield in the QGP scenario. In Fig.l, 
we have shown the ratio of photon yield obtained with 
one+two loop rate and one loop rate, for initial times Tj 
=0.2,0.4,0.6,0.8 and 1 fm. It can be seen for the lowest 
thermalisation time (tj=0.2 fm) two loop rate contribute 
more that 70% to the highp-r photons. If the higher loops 
contribute to the same order as the two-loop rate, their 
inclusion can change the two loop results drastically. In 
the given circumstances, it will be unwise to apply these 
results if the thermalisation time scale is as small as 0.2 
fm. 
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FIG. 1. Ratio of photon yield in QGP scenario obtained 
with one+two loop rate and one loop rate, for different initial 
times. The lines, from top to bottom corresponds to Ti=0.2, 
0.4, 0.6, 0.8 and 1 fm. Corresponding initial temperatures are 
listed in table 1. The lines are drawn with 

For Tj > .4 fm or more, the two loop rates contribute 
to 25% or less to the total yield. Thus for a realistic 
estimate of photon yield in the QGP scenario, one must 
consider Tj > .4 fm, when contribution of two loop rate is 
moderate and even if higher loops contribute to the same 
order, it will not affect the total yield substantially. Cal- 
culations with smaller Tj could not be relied upon. We 
thus assume Tj > 0.4fm. As will be shown later, for Tj > 
0.4fm, initially static QGP scenario could not explain 
the WA98 data. It under predicts them. However, if the 
fluid is assumed to have certain initial radial velocity, 
data could be fitted satisfactorily. Indeed, Peressounko 
et al H argued that the WA98 data could be fitted only 
with some initial radial velocity. Present analysis con- 
firms the result. Source of initial radial velocity could be 
the collisions among the constituents of the fluid (quarks 



and gluons). However, the high px photon yield will de- 
pend on the type of velocity distribution to be assumed. 
Peressounko et al J8| used linearly increasing velocity dis- 
tribution. In an earlier work || a Woods-Saxon form for 
the initial velocity distribution was used. Fig. 2, we have 
shown the computed photon yield for three types of ve- 
locity distributions, 
(a) Woods-Saxon type of velocity: 



v r (r) = u°/(l + exp({r - R„)/a) 



(b) Linearly increasing velocity: 
v r (r) -- 



v°r/R , r<R 
r > Rq 



(7) 



(8) 



and (c) surface peaked velocity: 

v r {r) = 4v° r exp((r - R )/a)/{l + exp{{r - R )/a)) 2 (9) 

with Rq =6.4 fm and a =0.54 fm corresponding to Pb 
nucleus. In these velocity profiles, the parameter wj? is 
chosen to be the maximum value of the distribution. We 
should remember that this choice of parametrization is 
not unique. For example, it may also be normalized to 
give the same average radial velocity. Due to the ambigu- 
ity of the definition of v®, the absolute values of photon 
yield for each profile in Fig. 2 may change in different 
normalization, but an important fact is that px depen- 
dence changes and the high pt part of the photon spectra 
is sensitive 
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FIG. 2. Photon yield for different types of initial velocity 
distributions, with t>°=0.2c. The dash-dot-dot line is for the 
surface peaked initial velocity, the dash-dot line is for the 
Woods-Saxon type of velocity distribution and the solid line 
is for the linearly increasing velocity profile. 

to the shape of the velocity distribution. Surface peaked 
velocity distribution is least effective for producing px- 



On the other hand, Saxon- Woods type of velocity pro- 
file produces more high px photons. Linearly increasing 
velocity lies in between. 

To see the physical picture of these profiles, let us con- 
sider some conserving quantity associated to the fluid. 
Then from the continuity equation for cylindrical sym- 
metric case, 

• 1 9 , ^ n 
n + (rnv r ) = 0, 

r or 

where n is the density of the conserved quantity. For 
example, if we take Wood-Saxon type velocity field, at 
r ~ 0, v r ~ const, so that 



-v r (0) 



dn 
Or 



near r ~ 0. Except for n (0) — > for r — > 0, the density 
variation at the center becomes infinitely large. Thus, 
the Wood-Saxon type profile represents a detonation at 
the origin, exploding from the center whole the material 
outwords. The linear profile corresponds to the homolo- 
gous expansion of the system, corresponding to a rather 
mild density change. For example, a typical solution is 



R 2 {ty \R(r) 



where / is an arbitrary function. On the other hand, 
when the QGP gas is initially almost at rest, and its 
surface area is expanding rapidly to the vacuum, similar 
to the early simple wave stage of the Landau model, then 
the velocity profile becomes surface peaked. If a nuclear 
collision generates initially almost static QGP fireball, 
the fluid velocity may be more surface peaked compared 
to the linear one. 

Before we compare our model prediction with the 
WA98 experiment, we would like to make few comments 
on the direct QCD and pre- equilibrium photons. Direct 
QCD photons are produced from the early hard collisions 
of partons in the nuclei and in Pb+Pb collisions make sig- 
nificant contribution to the high px yield. Gallmeister et 
al [0 claimed that prompt photons are able to explain 
the high px data in Pb+Pb collisions. Dumitru et al |2C| ] 
also arrived at a similar conclusion including the nuclear 
broadening effects. However, this point is still controver- 
sial due to uncertainties in prompt photon emission at 
AA collisions. Thus Alam et al H and also Srivastava 
and Sinha [^) calculated the prompt photon emission for 
Pb+Pb collisions. It was seen that for px >2 GeV, direct 
QCD photons alone can describe the data within a factor 
of 3-8 only. 

Pre-equilibrium photons are emitted before the estab- 
lishment of quasi equilibrium. Traxler and Thoma |2l| ] 
calculated pre-equilibrium photons and found them or- 
der of magnitude less than the equilibrium photons. Roy 
et al |E2] also calculated the pre-equilibrium photons. 



They used Fokker-Plank equations and found that pre- 
cquilibrium photons are less or at best equal to the equi- 
librium photons. The results refers to RHIC and LHC 
energies. At SPS one expects still lower contribution of 
pre-equilibrium photons. At SPS, the system is far from 
chemical equilibrium and number of quarks and anti- 
quarks will be less. Pre-equilibrium photons then can 
be neglected. 
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FIG. 3. Thermal plus hard QCD photons compared with 
WA98 experimental results. The initial time and temperature 
of QGP is 0.4 fm and 271 MeV respectively. The solid lines 
are obtained with initial radial velocity 0.2,0.3 and 0.4c (from 
bottom to top). The dotted line is obtained without any 
initial radial velocity. 

In Fig. 3 , we have compared the photon yield obtained 
in the present model with WA98 experiment. The initial 
time and temperature are Ti—OA fm and T,=271 MeV 
respectively. We have included the hard QCD photons 
for pr > 2 GeV following Alam et al B. Without any 
initial radial velocity (the dotted line) thermal photons 
together with hard QCD photons could not explain the 
data. High px part of the data are underestimated by 
a factor of 5 or more. It is possible to obtain a reason- 
able fit to data for still lower initial time(e.g. Ti=0.2 
fm, Ti=341 MeV), but as emphasized earlier, the results 
obtained thus are not reliable. Higher loops (three and 
more) contributions, which are neglected here may then 
be significant. When we suppose that the initial tem- 
perature is not so high, then the data requires an initial 
radial velocity in our model. In Fig. 3, the solid lines 
are the photon yield obtained with initial radial velocity 
with maximum velocity i>J?=.2,.3,.4c. Here we take the 
surface peaked initial velocity profile, which corresponds 
to a scenario of QGP fireball with surface expansion. As 
expected with initial v r , high px photon yield is increased 
and we find that the data is well describe if v° lies in 
the range of 0.3-0. 4c. When we assume more mild QGP 



gas expansion represented by the linear velocity profile, 
then these values are little bit smaller. In Fig. 4, same 
results are shown for initial time t,;=1.0 fm. In this case, 
without any initial radial velocity, data are more under- 
predicted. Reasonable agreement with data is obtained 
again with (surface peaked) initial velocity in the range 
of vj?=0.4-0.5c. The present results are in agreement with 
the finding of Peressounko et al || that the WA98 single 
photon data require initial radial velocity of the fluid, at 
least in the surface region. 
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FIG. 4. Thermal plus hard QCD photons compared with 
WA98 experimental results. The initial time and temperature 
of QGP is 1.0 fm and 200 MeV respectively. The solid lines 
are obtained with initial radial velocity 0.3, 0.4 and 0.5c (from 
bottom to top). The dotted line is obtained without any 
initial radial velocity. 

We now show the prediction for RHIC energy. In Fig. 5, 
we have shown the results obtained for Au+Au collisions 
at y / s=130 GeV. The initial time and temperature was 
assumed to be r,=0.5 fm and Ti=300 MeV ||. The 
dash-dot line is the hard QCD photons at RHIC energy. 
Pt broadening effect is ignored. At RHIC energy the ef- 
fect is not as important as at SPS energy. The dotted 
lines are the thermal photons obtained for initial radial 
velocity v®=0,0.3 and 0.5c. Solid lines show the total con- 
tribution ( thermal and hard QCD photons). It can be 
seen that for pr < 2 GeV, thermal photons dominate the 
spectrum. Hard QCD photons dominate at higher pt > 4 
GeV. In the intermediate pr range (2 GeV-4 GeV) , ther- 
mal photons contribute as copiously as the hard QCD 
photons. This is the range where initial fluid velocity 
has considerable effect. Considering the importance of 
hard QCD photons at RHIC energy, it is of utmost im- 
portance to determine experimentally the single photon 
spectra in pp collisions at RHIC energy. Only then it will 
be possible to comment on the possible QGP formation 
at RHIC energy from the single photon data. 



To summarize, we have analyzed the recent WA98 sin- 
gle photon data assuming QGP formation in the initial 
state. It was shown that in order to have a reliable es- 
timate of photon yield, initial time of the hydrodynamic 
evolution should be larger than 0.4 fm, otherwise neglect 
of three and higher loops in the photon production rate 
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FIG. 5. Theoretical estimate of direct photons at RHIC 
energy for Au+Au collisions at B=130 GeV. The dash-dot line 
is the hard QCD photons. The dotted lines are the thermal 
photons with initial radial velocity uj?=0,0.3 and 0.5c. The 
solid lines are the total contribution of hard QCD photons 
and thermal photons. 

would be unjustified. They can contribute in the same 
order as the two loop photon rate (which is more than 
70% at high pr for an initial time of 0.2 fm). As far as 
Tj > 0.4 fm, then initially static QGP scenario underes- 
timates the high pt photons compared to the observed 
data. The data could be fitted only if some kind of ini- 
tial radial velocity is present. Photon yield depend on the 
initial velocity distribution, changing by factor of 2-3 at 
the high px , from the surface emission scenario (surface 
peaked) to the detonation at the center (Woods-Saxon 
form). WA98 data could be fitted reasonably well if the 
thermalisation time Tj, ranges between 0.4-1.0 fm (corre- 
sponding to initial temperature between 271-200 MeV), 
and if the QGP fluid has an initial (surface peaked) radial 
velocity in the ranges of 0.3-0. 5c. The model was applied 
to predict photon spectra at RHIC energy in Au+Au col- 
lisions. It was shown that high px spectra is dominated 
by hard QCD photons. Thermal photons dominate low 
Pt part of the spectra and in the intermediate range of pt 
thermal photons contribute comparably to hard photons. 
We feel it is important to measure hard QCD photons in 
pp collisions at RHIC energy. Only then it will then be 
possible to comment on possible QGP formation from 
the direct photon data at RHIC. 
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TABLE I. The initial temperature of the QGP for certain 
initial times n. dN/dY was assumed to be 750. 



Ti(fm) 


0.2 


0.4 


0.6 


0.8 


1.0 


T? GP (MeV) 


341 


271 


237 


215 


200 



